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Why are central radio relics so rare? 



o 

> 
O 



O 

u 

o 



> 

o 



X 



F. Vazza^'^*, M. Briiggen^, R. van Weeren^, A. Bonafede^, K. Dolag*, G. Brunetti^ 

^ Jacobs University Bremen, Campus Ring 1, 28759, Bremen, Germany 

^ INAF/Istituto di Radioastronomia, via Gohetti 101, 1-40129 Bologna, Italy 

^ Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands 

^ University Observatory Munich, Scheinerstr. 1, D-81679 Munich, Germany 



Accepted ???. Received ???; in original form ' 



ABSTRACT 

In this paper we address the question why cluster radio relics that are connected to shock 
acceleration, so-called radio gischt, have preferentially been found in the outskirts of galaxy 
clusters. By identifying merger shock waves in cosmological grid simulations, we explore 
several prescriptions for relating the energy dissipated in shocks to the energy emitted in the 
radio band. None of the investigated models produce detectable radio relics within 100-200 
kpc from the cluster centre. All models cause > 50 per cent of the detectable relic emission 
at projected distances > 800 kpc. Central radio relics caused by shocks that propagate along 
the line-of-sight are rare events for simple geometrical reasons, and they have a low surface 
brightness making them elusive for current instruments. Our simulations show that the radial 
distribution of observed relics can be explained by the radial trend of dissipated kinetic energy 
in shocks, that increases with distance from the cluster centre up until half of the virial radius. 
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1 INTRODUCTION 



Radio relics can be divided into two main groups jKempner et al.l 
I2OO4I) : Radio gischt are large elongated, often Mpc-sized, radio 
sources located in the periphery of merging galaxy clusters. In 
a few cases, their spatial co-location with shocks in the ther- 
mal gas has been dete rmined (e.g. in the case of Abell 3667, 
iFinoguenov et alj|20ict) . They most likely trace shock fronts in 
which particles are accelerated via diffusive shock acceleration 
(e.g. Blandford & Eichler 1987). Among them are double-relics 
with t he two relics located on different sides of a cluster centre 
(e.g., iBonafede et al . 'v an Weeren et a i] |2009l : IVenturi et al.l 



2007llBagchi et al 



2006; Rottg ering et al.ll997l : lBrown et alj20l'll : 
Bagchi et alj 1201 1 ). Radio phoenices are related to radio-loud 
AGN. Fossil radio plasma from a previous episode of AGN activity 
is thought to be compressed by a merger shock wave which boosts, 
both, the magnetic field inside the plasma as well as the momenta 
of the relativistic particles. Here we are only concerned with radio 
relics of the gischt type. 

The sizes of relics and their separations from the cluster cen- 
tre vary significantly. Examples for radio relics wi th sizes of 1 Mpc 
or ev en larger have been observed in Coma jGiovannini et al.l 



199 ih. Abell 2255 Jperetti et alj Il997h and Abell 2256 (see 



van Weeren et al.f201 iL for other examples and references). 



Even though the sample of known radio relics is still small 
and incomplete (34, see Table 1), one can start to find correla- 



tions between size, location and spectral index of these unique 
sources, which can be compared to simulations of relic forma- 
tion, van Weere n et al. j 12009) find that on average the smaller 
relics have steeper spectra. Such a trend is in line with pre- 
dictions from shock statis t ics derived from cosmological simu- 
lations jVazza et alj |2()10| ; Iskillman etd] I2OO8I : IPfrommer et al.l 
l2007l : lHoeft et al.ll2008 l).Thev find that larger shock waves occur 
mainly in lower-density regions and have larger Mach numbers, 
and consequently flatter spectra. Conversely, smaller shock waves 
are more likely to be found in cluster centres and have lower Mach 
numbers and steeper spectra. In this paper we explore why radio 
gischt relics are never observed close to the centre of galaxy clus- 
ters, using a set of high-resolution cosmological re-simulations of 
massive galaxy clusters at three cosmic epochs (z=0.0, z=0.3 and 
z=0.6). 



2 SIMULATIONS OF RADIO EMISSION FROM SHOCKS 

The simulations analyzed in this work were produced with the 
Adaptive Mesh Refinement code ENZO 1.5, developed by the Lab- 
oratory for Computational Astrophysics at the University of Cal- 
ifornia in San Diego 0. A detailed discussion on the numerical 
se tup adopted to produce the simulations of this work is presented 
in IVazza et alj JIoiQ) . Twenty galaxy clusters with masses in the 
range 6 X 10" s; M/Mq < 3 X lO^'^ were extracted from a total 
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Figure 1. Emitted "radio" power, P, from simulated relics (colours, in units of logjg P/ [erg/s]), for the same galaxy cluster at z=0 and for 3 emission models 
(KJ, KJ JiparallelRadial and KJJitangledRadial, in the last two models we assumed Bq = 5 /iG). The side of each image is 10 Mpc/h. 



cosmic volume of Lbox ~ 480 Mpc/h. With the use of a nested 
grid approach to produce initial conditions at high resolution in the 
region of cluster formation, the final mass resolution for the Dark 
Matter (DM) component is mdm ~ 6.76 ■ 10^ Mq. The refinement 
criterion is based on gas/DM over-density as well as ID velocity 
jumps which results in a peak resolution of 25 kpc/h inside a radius 
of ~ 5 — 6 Rv (Ry begin the virial radius) from the centre of each 
cluster (see lVazza et al.l l2010' for further details). 

The assumed cosmology is the "concordance" ACDM one, 
with Qo = 1-0, = 0.0441, ^dm = 0.2139, Qa = 0.742, 
Hubble parameter h — 0.72 and a normalization for the primordial 
density power spectrum erg = 0.8. Our simulations neglect radia- 
tive cooling, star formation and AGN feedback processes. In this 
work we analyse cluster s at redshifts z=0.0, z=0.3 and z=0.6 (see 
IVazza et al .I2OIOI I20T Tal for complementary studies of our clusters 
at these epochs). 

2.1 Shock detection 

The shock waves in our r uns are detected in post-processing with 
the algorithm presented in lVazza et iDiooi): this scheme is based 
on the analysis of ID velocity jumps across cells. The 3D diver- 
gence of the vel ocity marks the centre of the shock region (see also 
IRvu et al.ll2003h . The one-dimensional Mach number for each cell 
boundary is reconstructed by inverting the ID velocity discontinu- 
ities along each scan axis using the Rankine-Hugoniot jump rela- 
tions. The final Mach number is reconstructed by combining the 
three ID jumps. This method has been extensively tested ag ainst 
similar methods used in grid codes (e.g. IVazza et al.|[2011bl) . and 
has proven to be an effici ent and accurate m easure of shock waves 
in cosmological runs (see l Vazza et al.l2010l . for the statistical prop- 
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Figure 2. Sample of the 6 brightest "relics" in our sample 
{KJ JiparallelRadial model assumed, with Bo = 5 ^iG). The colours 
represent the projected X-ray bolometric emission (in arbitrary counts), the 
contours represent the objects identified with our relic finder. The size of 
each image is 3.5 X 2.2 Mpc/h. 
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Table 1. List of gischt radio relics. From left to right: name of host cluster, 
redshift and projected distance from the centre. In case of multiple relics per 
cluster, the letters N,S,E,W denote the direction with respect to the cluster 
centre. We notice that in the case of a few objects (labelled with "*" in the 
table) the identification as "radio gischt" can be debated. 



name 


z 


r[kpc] 


Coma cluster 


0.023 


2008 


Abell 548b 


0.042 


693 


Abell 548b* 


0.042 


456 


Abell 3376 E 


0.0456 


845 


Abell 3376 W 


0.0456 


1008 


Abell 3667 E 


0.055 


1459 


Abell 3667 W 


0.055 


1887 


Abell 2256 


0.059 


470 


CIZA J0649.3+1801 


0.064 


802 


RXC J1053.7+5452 


0.07 


981 


Abell 2061 


0.078 


1558 


Abell 2255 


0.0806 


1065 


Abell 3365 E 


0.093 


1079 


Abell 3365 W 


0.093 


777 


ZwCl 0008.8+5215 E 


0.104 


798 


ZwCl 0008.8+5215 W 


0.104 


808 


Abell 1664* 


0.128 


895 


Abell 1240 S 


0.159 


982 


Abell 1240 N 


0.159 


924 


Abell 2345 E 


0.177 


1431 


Abell 2345 W 


0.177 


857 


Abell 1612 


0.179 


894 


C1ZAJ2242.8+5301 N 


0.192 


1565 


CIZAJ2242.8+5301 S 


0.192 


1062 


Abell 115 


0.197 


998 


Abell 2163 


0.203 


1405 


Abell 746 


0.232 


1606 


RXCJ1314.4-2515 W 


0.244 


577 


RXC J13 14.4-25 15 E 


0.244 


937 


Abell 521 


0.247 


692 


ZwCL 2341.1+0000 N* 


0.27 


1198 


ZwCL 2341.1+0000 S 


0.27 


767 


Abell 2744 


0.308 


1501 


MACSJ0717.5+3745* 


0.555 


298 



erties of shocks in the same sample of clusters used in this work). 
The co-moving kinetic power through each shock surface is given 
by 

Fke = ^AS, (1) 

where pu is the co-moving up-stream density, Vs = Mcs is the co- 
moving speed of the shock, M is the upstream Mach number and 
AS is the surface area of the shocked cell. 



2.2 Radio emission from shocks 

For each cluster, we generated mock radio observations along each 
of the tftree coordinate axes by creating maps of Fke, weighted by 
a function, /m. This function encapsulates the poorly constrained 
dependence between the Mach number of the shock and the re- 
sulting radio emission. For instance, the acceleration efficiency of 
cosmic rays protons at shocks, the electron-to-proton energy ratio, 
the amplification of the magnetic field due to CR instabilities are 
still largely unknown (e.g.lBri iggen et al.l201ll . for a recent review). 
Here we explore 6 different forms of /m . 

The model Kin uses /m = 1, and illustrates the simple trend 



of the kinetic power at shocks in the simulated volume (in this case, 
the acceleration efficiency of particles is independent of M). 

In the model KJ w e use t he Mach number-dependent effi- 
ciency of acceleration o f iKang & Jones (2007), /m = niM), fol- 
lowing the results of ID studies of Diffusive Shock Acceleration 
jKang & Jonesll2007l) Pl 

In a second set of models, we include the radial decrease 
of the cluster magnetic field and the amplification of the pre- 
shock magnetic field by shocks. The radial model of reference 
for the "background" clu ster magnetic field is the one inferred 
by iBonafede et al.l feOlOl) for observations of the Coma cluster: 
B{p) = i?o(p/p250o)", with a « 0.5 and P2500 the density of 
the cluster core. 

In the model KJJitangledRadial, we include a factor, R — 
Bd/Bu = A'P /{M^ + 3), for the amplification of an isotropic 
upstream magnetic field following the compression by the shock. 

At cluster shocks, however, the amplified magnetic fields may 
not be isotropic and thus we also explore an alternative model. A 
magnetic flux tube of strength B^, oriented at an angle 6^ with 
respect to the shock normal will be bent and amplified, according 
to 

Bd = Bu (cos^ 6»u + sin' 6^) ^''^ , (2) 

where the suffices u and d refer to the up- and the downstream re- 
gions, respectively, and i? = pd/pu is the shock compression ratio. 
From equation ([2} it is easy to see that the magnetic field amplifica- 
tion is negligible when B^ is oriented parallel to the shock normal 
(6'u ~ 0, thereby B^ — Bu), whereas B^ — RBu when Bu is 
nearly parallel to the shock surface. In models KJ JiparallelRadial, 
we thus assume that the amplified magnetic field at the shock is 
Bd ~ RBu- 

In all cases the final weighting function for the kinetic en- 
ergy flux is /m = rj{M){RB{p)f/[{RB{p)f + Bf^], where 
Bic ~ 3.2 ■ (1 + z)^fiG is the magnetic field equivalent to 
the IC losses from CMB photons, at each redshift). In the model 
K JJitangledRadial we a dopted the central value of Bo = 5 pG as 
m iBonafede etai] ilOld) . while in the models KJ JiparallelRadial 
we used Bo = 1 pG, = 5 pG and = 10 pG. We assume here equal 
average values of the magnetic field for relaxed and very perturbed 
objec ts, which is supported by recent observations IBonafede et^ 
l20Ilh . The formulae for /m bracket reasonable physical models, 
but they cannot provide exact absolute values of the radio power in 
the relics. In Sec l2.3l all radio maps are normalized to the maximum 
within the image. 

In each shocked cell of our 3D grid, we compute the energy 
that is thus converted to radio emission using the locally computed 
values of M, p and Vs- This emission is then projected along each 
of the three coordinate axes, for a total of 60 mock radio images for 
each assumed model and at each r edshift. Generatin g projections 
along arbitrary angles of view (e.g. lHoeftet ^|2008|) is computa- 
tionally non-trivial because the computational boxes contain a large 
number of cells (> 500"^), and, hence, we defer it to future work. 
The pixels in our images have the same size ( comoving 25 kpc/h) 
of the maximum available 3D resolution, corresponding to the size 



^ The details of particle acceleration in the regime of Mach numbers typ- 
ical of the ICM, 2 ^ M ^ 5 are not yet well-constrained. For instance, 
recent studies with particle-in-cells methods investigated additional accel- 
eration mechanism for particles at shocks (e.g. shock drift acceleration), 
suggesting the possibilities also of a different trend with Mach number, an d 
large efficiency for low Mach number shocks jOargate & Spitkovskvl201ll) . 
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of the primary beam of the Very Large Array in D configuration 
(45" X 45") for a luminosity distance of « 220 Mpc (z^ 0.05). 

In Figure[T]we show mock radio emission maps for one simu- 
lated cluster at z=0, for 3 investigated emission models. This clus- 
ter experienced a major merger ~ 0.5 Gyr ago, and two powerful 
merger shocks are expanding out of the main cluster core. These 
shocks are ~ 1 Mpc large and morphologically similar to real ob- 
served relics, and they are clearly visible in all investigated emis- 
sion models. The much softer envelope of internal merger shocks 
and of some large-scale expanding outer ones is very difficult to de- 
tect with current radio telescopes, given that their power is several 
orders of magnitude dimmer than the brightest double relics, for all 
investigated models (see SectO. 

2.3 Relic detection 

The mock radio observations at all redshifts were analysed with a 
2D relic detection scheme, which identifies as part of the same relic 
all pixels brighter than a given threshold. Here we do not aim to re- 
produce the spectral and absolute power properties of relics, but 
only focus on their geometrical properties in projection. We thus 
normalize the pixel values of each map to the maximum within the 
image, and we run our relic finder over the 10 per cent most lu- 
minous pixels in each map. In most cases, the brightest pixels are 
connected to form relatively large (> 200 comoving kpc) objects. 
However, depending on the threshold, some of the them can be dis- 
connected by a few cells from the main structure. We found that by 
considering all pixels detached by less than 8 pixels in each direc- 
tion as part of the same structure, we obtain a reasonable identifica- 
tion of what radio observers would typically consider a relic (tests 
with smaller linking lengths between pixels, 2 and 4, have yielded 
nearly identical results). Once we have identified all radio relics in 
each image and for every emission model, we determine their total 
power, size and distance from the centre of the host cluster (com- 
puted as the projected distance from the radio luminosity weighted 
centre of each relic and the X-ray peak). 

In Figure |2l we show the contours of the six most pow- 
erful relics detected in our sample of 60 maps at z=0 for 
KJ JiparallelRadial with Bo = 5 /iG, as detected by our algorithm. 
The first three panels show the emission from the same object as 
in Fig[T]seen along the three coordinate axes. In the other panels 
the emission comes from different simulated galaxy clusters. In all 
panels, we also over-plotted of the X-ray brightness map for each 
cluster in colours. 



3 RESULTS 

The main result of this paper is presented in Fig[3] Here we show 
the cumulative distributions of projected comoving distances from 
the centre for the relics from our sample (Table 1) as well as the 
radio relics detected in our mock images for 6 different emission 
models at the three epochs. Since we do not intend to reproduce 
the absolute power of relics (and their chance of detection in mock 
images), we normalize all distributions to the observed one at a 
comoving radius of 2 Mpc. 

In the observed distribution, ~ 50 per cent of relics are located 
at distances > 800 kpc from the cluster centre. In only one object 
of the sample (~ 3 per cent), the projected radial distance is below 
300 kpc (MACSJ07 17.5+37450). At all epochs our simulated dis- 



tributions show a very similar trend, with the exception of model 
KJ, in which the distributions are more extended 50 per cent of 
the relics is found at r > 1000 kpc). Notably, including a profile 
for the magnetic field distribution shifts the distribution to slightly 
lower radii compared to the other models. The resulting distribu- 
tion matches the observed distribution in the range r ^ 500 — 2000 
kpc quite well. Comparing the simulated and observed distribution 
at larger redshifts is more difficult, since only a few observations 
are available beyond z~0.3. However, Fig|3] shows that the radial 
distribution of simulated relics at higher redshifts is more radially 
concentrated compared to z=0, due to the fact that clusters are phys- 
ically smaller, and that major merger are more frequent at large 
redshifts. In all cases, even at higher redshifts the frequency of cen- 
tral radio relics is very low, ^ 5 percent for r 300 kpc for all 
investigated models. 

The observed distribution of relics for r > 800 kpc at z=0 
is best reproduced by models that include the spatial variation of 
the background cluster magnetic fields, and the local amplification 
by shocks. Interestingly, though, the cumulative distribution is al- 
ready apparent in a model where /m = 1. This implies that the 
radial distribution of radio gischt relics is driven by the radial dis- 
tribution of Fke. Hence, we have plotted radial profiles of Fke, 
volume-averaged within radial shells, for all simulated clusters at 
z=0 (top panel of Fig|4]l. Fke initially increases with radius, be- 
cause the surface of the relics, AS, increases roughly as r^, and 
Vg also increases with radius because of the radial drop of the ICM 
temperature outwards. Eventually though, the decrease in pu (as 
~ for a /3-model outside of the cluster core, /? — 2/3) wins 
over the increase in AS, and the profile of Fke will turn over. The 
maximum of the distribution, however, depends on the assumed 
/m- For a fixed acceleration efficiency, the peak is at ^ 400 — 500 
kpc, and increases to up to 1 Mpc if the dependence of the ac- 
celeration efficiency on Mach number, ri{M), is included. In all 
models, we expect the number of radio relics to go down again 
with cluster-centric distance beyond a radius of ~ 0.5 — 1 Mpc. 
Since our sample has a mass range of 0.6 dex, we also computed 
the profiles normalized to the respective virial radii of the cluster 
(bottom panel of Fig|4j, and found a very similar behaviour, with 
maxima at ~ 0.2 — 0.5J?v. 

To test the role played by geometrical effects, we over-plot 
in the same figure the result of very simple 2D geometrical model 
(^'geometrical model 1 "), which assumes that all radio relics form 
in a major merger with the initial size of the core, rcore (we fix 
Tcoro ~ 350 has an average value for our clusters). The radio relics 
increase their linear size while expanding the outer cluster atmo- 
sphere, as Zreiic('') « r — Tcoro. The relics are assumed to be visible 
only if the plane of the merger is in the plane of the sky, and are as- 
sumed to be invisible elsewhere due to the lack of limb brightening. 
In this case, the cumulative probability of finding a radio relic at a 
given distance from the main cluster centre is given by the ratio be- 
tween the size of the relic at a given radius and the circumference of 
the corresponding radial shell, A'^(> r) = Z(r)/r = (r — rcoic)/r. 
We find that this geometric model gives a reasonable description 
of the observed relic distribution for r > 800 kpc (thus suggest- 
ing that the effect of limb brightening is very important for the 
outermost relics), but that it offers only a poor description of the 
observed relic distribution at smaller radii. 

One may wonder why one never observes as relics the shocks 



^ As noted in Tab. 1 , the identification as a "radio gischt" is also debatable. 



and the location of the X-ray centre of the host cluster is uncertain (see 
iBonafede et al.ll2009l;lvan Weeren et alj2009l) . 
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Figure 3. Cumulative distributions of simulated radio gischt emissions for the different emission models at z=0, z=0.3 and z=0.6. The dashed dark grey line is 
for observed radio relics; the light grey lines are for the geometrical models 1 and 2 

, outlined in the text (see Sec[3]for details). 



in the mid-plane that foim when clusters collide before core pas- 
sage. These shocks turn out to have very low Mach numbers be- 
cause a filament connects the two clusters. The diffuse gas in these 
filaments is already > lO*" K and this gas gets adiabatically com- 
pressed and hence heated as the two clusters approach each other. 
Hence, we do not find these shocks as radio relics in our cosmolog- 
ical simulation. 

We note that the probability of observing radio relics in projec- 
tion against the cluster centre is very low for geometrical reasons. 
Even assuming that the occurrence and properties of radio relics 
are independent of distance, only a small fraction of relics should 
appear projected onto the cluster centres. This follo ws from simple 
geom etrical considerations (see Suppl. Material in ISrunetti et al.l 
120081 "geometrical model 2" in Fig[3}: relics projected to within 
distance a r from the cluster centre are those that lie along a line- 
of.sight in a cylinder of radius r, that intercepts a decreasing frac- 
tion of the cluster volume with radial distance. The resulting distri- 
bution of probability is thus: 

F{< r) = / sin(6l)d6l / I'^dl, (3) 

^2Mpc Jo Jo 

where a{6) — min(r/ sin(6'), i?2Mpc), i?2Mpc — 2 Mpc. 

The resulting distribution of relics is shown in Fig|3l the dim- 
ming of projected relics due to the lack of limb-brightening will 
further decrease the chance of observing them at small distances. 

From the comparison with our simulated relics, it can be seen 
that the effect of dimming further kills the probability of detecting 
relics in projections on the cluster core, even in those few cases ex- 
pected by geometry. Indeed, even in the case of the most powerful 
relic (Fig[T) merely a bright spot of emission is detected at ~ 100 
kpc from the centre in the j/-projection, but most of the extension of 
the two giant relics is lost because of the lack of limb brightening 
in this case. The flux goes down by the ratio between the largest 
linear size of relics to their widths, which is typically > 10, (e.g. 
iHoeft et alJl2008l : lBums & Skillmanl2oTll) . 

Based on our maps, it is likely that relics in projection against 
the centre of the cluster have not been detected so far because of 
sensitivity limits of radio observations. When a shock is seen face- 
on, and the resulting relic is projected against the cluster centre, the 
peak of its radio emission lies at the edge of the cluster; its the radio 
brightness would be ^ 100 times lower in projection onto the clus- 
ter centre. Relics have a low surface brightness (~ ifJ-Jy arcesc"^ 



at 1.4 GHz) and apart from a few cases, are detected at a few a 
by current radio telescopes. Hence, detecting central relics as those 
seen in our runs would require a sensitivity currently unachievable, 
also for future radio telescopes (such as LOFAR and SKA) because 
of confusion limits. 



4 CONCLUSIONS 

In this paper we explored the question why gischt radio relics are 
never detected in the innermost regions of clusters. Using simple 
prescriptions to model the synchrotron emission emitted at shock 
waves, we produced mock radio maps of a sample of 20 massive 
galaxy clusters simulated at high resolution using the ENZO code, 
at redshifts z=0.0, z=0.3 and z=0.6. Our models parametrized the 
efficiency with which electrons are accelerated at shocks, the back- 
ground magnetic field and the amplification of the upstream mag- 
netic field. In all models, no detectable radio-emission from radio 
relics is found inside the innermost ~ 100 — 200 kpc of our clus- 
ters. Nearly ~ 50 per cent of detectable radio relics are located at 
projected distances > 800 kpc (> 0.2 — 0.5i?v) from the cluster 
centres. The low probability of observing radio relics close to clus- 
ter centres results from the fact that the kinetic energy dissipated in 
merger shock waves increases with distance from the cluster cen- 
tre, combined with the more obvious geometrical reasons. This is 
because the Mach number of merger shocks as well as their surface 
areas increase with distance from the centre faster than the density 
decreases. The acceleration efficiency as a function of Mach num- 
ber as well as the radial dependence of the cluster magnetic field do 
not alter this trend qualitatively. Finally, the lack of limb brighten- 
ing makes the detection of the few central projected relics almost 
impossible at the typical sensitivity of all existing radio telescopes. 

In summary, the lack of central radio relics in observations 
is due to the low energy dissipation of shock waves in cluster cen- 
tres. Attempts to constrain the mechanism at work in the production 
of radio reli cs have been already explored in cosmological simu- 
lations (e.g.lMiniati et al.ll200ll ; faoeift et alJboOSi ; IPfrommer etaD 
|2007| ; ISkillman et al.ll201 ll) . In the future we also plan to study in 
more detail the statistical, spectral and morphological properties of 
simulated radio relics in our sample. 
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Figure 4. Volume-averaged (spherically and over all clusters at z=0) pro- 
files of the kinetic power at shocks, Fke, for the emission models inves- 
tigated (for clarity only model KJ JipamllelRadial with Bq = 5 /xG is 
shown). Top panel: profiles averaged at physical radii. Bottom panel: Same 
profiles as in the top panel, only as a function of radius normalized to the 
virial radius, iJv, of each cluster. 
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